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PREFACE 


The  investigation  reported  herein  was  authorized  by  the  Long  Beach 
Port  Authority  under  a contract.  Agreement  No.  WES  76-U,  dated  12  Sep- 
tember 1975. 

The  investigation  was  conducted  from  September  1975  to  October  1976 
by  personnel  of  the  Hydraulics  Laboratory,  U.  S.  Army  Engineer  Waterways 
Experiment  Station  (WES),  under  the  direction  of  Mr.  H.  B.  Simmons, 

Chief  of  the  Hydraulics  Laboratory,  Dr.  R.  W.  Whalin,  Chief  of  the  Wave 
Dynamics  Division,  and  Mr.  C.  E.  Chatham,  Chief  of  the  Harbor  Wave 
Action  Branch.  Dr.  J.  J.  Wanstrath,  Research  Oceanographer,  conducted 
the  investigation  and  prepared  this  report.  Mr.  R.  R.  Bottin,  Jr., 
aided  in  the  development  of  the  finite  element  grids.  Drs.  H.  S.  Chen 
and  C.  C.  Mei  of  the  Massachusetts  Institute  of  Technology  provided 
documentation  of  the  hybrid  finite  element  computer  program  they  de- 
veloped and  materials  to  aid  in  its  utilization. 

Directors  of  WES  during  the  investigation  and  the  preparation  and 
publication  of  this  report  were  COL  G.  H.  Hilt,  CE,  and  COL  John  L. 
Cannon,  CE.  Technical  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  converted 
to  metric  (Si)  units  as  follows: 


Multipl; 


To  Obtain 


0.301+8 

1+01+6.856 

0.761+551+9 

0.301+8 

0.09290301+ 


square  metres 
cubic  metres 


cubic  yards 
feet  per  second 
square  feet  per  second 


LONG  BEACH  HARBOR  NUMERICAL  ANALYSIS  OF  HARBOR  OSCILLATIONS 


ALTERNATE  PLANS  FOR  PIER  J COMPLETION  AND  TANKER 


TERMINAL  PROJECT  (NO  LANDFILL) 


PART  I : INTRODUCTION 


Objective 


1.  The  objective  of  this  study  was  to  investigate  by  use  of  a 
numerical  model*  the  response  to  long-period  wave  excitation  of  a pro- 
posed modification  to  Pier  J in  Long  Beach  Harbor  for  oil  tanker  berth- 
ing and  general  cargo  facilities.  Wave-height  amplification  factors  ar 
normalized  maximum  current  velocities  were  plotted  versus  wave  period  a 
the  three  proposed  tanker  terminals  to  ascertain  whether  or  not  signifi 
cant  harbor  oscillations  might  occur  that  would  pose  problems  to  ship 
berthing. 


Scope 


2.  The  tanker  terminal  configuration  under  consideration  in 
this  study  (plan  3)  omits  the  landfill  of  approximately  111  acres** 

(11  million  cu  yd)  on  the  end  of  Pier  J (Figure  l).  Three  other  plans 
have  been  testedt  and  are  summarized  on  the  following  page  (Figures  2 


* H.  S.  Chen  and  C.  C.  Mei,  "Oscillations  and  Wave  Forces  in  an  Off- 
shore Harbor  (Applications  of  the  Hybrid  Finite  Element  Method  to 
Water-Wave  Scattering),"  Report  No.  190,  197*+,  Massachusetts  Institu 
of  Technology,  Cambridge,  Mass. 

**  A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (Si)  units  is  presented  on  page  3* 
t J.  R.  Houston,  "Long  Beach  Harbor  Numerical  Analysis  of  Harbor 
Oscillations;  Alternate  Plans  for  Pier  J Completion  and  Tanker  Termi 
nal  Project,"  Miscellaneous  Paper  H-76-20,  Report  2,  Sep  1976, 

U.  S.  Army  Engineer  Waterways  Experiment  Station,  CE,  Vicksburg, 
Miss. 


Plan 


1A 


IB 


2 


3 


Designation 


Description 


Scheme  STFP  2 


Scheme  STFP  2 


Scheme  STFP  3 


300- ft  opening  (1000-ft-long  east 
breakwater)  at  the  east  end  of  the 
62-ft-deep  basin  requiring  111  acres 
of  landfill  on  Pier  J. 

600-ft  opening  (700-ft-long  east 
breakwater)  at  east  end  of  the 
62-ft-deep  basin  requiring  111  acres 
of  landfill  on  Pier  J. 

1100- ft  opening  (1300-ft-long  east 
breakwater)  at  the  east  end  of 
the  62-ft-deep  basin  requiring 
50  acres  of  landfill  on  Pier  J. 

1700-ft  opening  (700-ft-long 

east  breakwater)  at  the  east  end 
of  the  62-ft-deep  basin  requiring  no 
landfill  on  Pier  J. 


3.  For  all  plans,  including  plan  3,  the  incident  waves  were  di- 
rected normal  to  the  main  breakwater  with  periods  from  30  sec  to  300- 
360  sec.  For  this  study,  wave  amplitudes  and  currents  were  calculated 
every  6 sec  for  the  period  range.  Resonant  peaks  also  were  defined  by 
considering  incident  wave  periods  in  increments  as  low  as  1 sec. 
Frequency  response  curves  are  given  comparing  plans  IB  and  3 in  detail, 
and  a tabular  summary  of  resonant  peaks  is  presented  for  all  plans. 


) 
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LANDFILL  AREA 


TERMINAL 


TERMINAL 


TERMINAL 


BREAKWATER 


Location  of  tanker  terminals 


Figure  2.  Scheme  STFP  2 


PART  II:  NUMERICAL  MODEL 


. The  response  of  the  tanker  terminal  basin  to  long  wave  excita- 
tion was  determined  by  using  a hybrid  finite  element  numerical  model 
developed  recently  by  Chen  and  Mei.*  The  model  solves  the  following 
generalized  Helmholtz  equation: 


V • [h(x,y)V<j>(x,y)]  + — <j>(x,y)  = 0 (l) 

where  <t>(x,y)**  is  the  velocity  potential  defined  by  u(x,y)  = -V<f>(x,y), 
with  u(x,y)  being  a two-dimensional  velocity  vector  and  w an  angular 
frequency.  Equation  1 governs  small  amplitude  undamped  oscillations 
of  water  in  a basin  of  arbitrary  shape  and  variable  depth  forced  by 
periodic  long  waves.  It  has  been  further  assumed  that  the  flow  is 
irrotational . 

5.  The  boundary  condition  along  the  shoreline  and  along  the  de- 
tached breakwater  surrounding  the  three  proposed  berths  is  that  the 
normal  component  of  the  velocity  be  equal  to  zero.  Therefore,  the 
breakwater  is  considered  as  a solid  barrier.  No  special  boundary  con- 
ditions are  made  to  account  for  any  vessels  utilizing  the  tanker  termi- 
nal basin,  since  for  long-period  waves  of  small  amplitude  ship  heave 
equals  incident  wave  height.  Therefore,  the  vessels  do  not  signifi- 
cantly alter  the  oscillation  characteristics  of  the  basin. 

6.  The  Helmholtz  equation: 

? 2 

V <j>(x,y)  + ~r  <t>(x,y)  = 0 (2) 


is  the  governing  equation  for  a constant-depth  ocean  region  outside  the 


basin. 


7 • For  a harbor  in  a semi-infinite  ocean  with  a straight  coastline 


* Chen  and  Mei,  op.  cit. 

**  For  convenience,  symbols  and  unusual  abbreviations  are  listed  and 
defined  in  the  Notation  (Appendix  A). 
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there  is  an  incident,  reflected,  and  scattered  wave.  The  scattered 
wave  has  a velocity  potential  <t>s  given  by 


d>  = Y a H (kr)  cos  n9 
s n n 

n=0 


where  a are  unknown  coefficients  and  H (kr)  are  Hankel  functions 
n n 

of  the  first  kind  of  order  n . 4>  satisfies  the  radiation  condition 

s 

that  the  scattered  wave  must  behave  as  an  outgoing  wave  at  infinity. 
This  condition  is  known  as  the  Sommerfeld  radiation  condition  and  may 
be  expressed  mathematically  as  follows: 


~ ik)  = 0 


8.  Chen  and  Mei  used  a calculus  of  variations  approach  and  ob- 
tained a Euler-Lagrange  formulation  of  the  boundary  value  problem.  The 
following  functional  with  the  property  that  it  is  stationary  with  re- 
spect to  arbitrary  first  variations  of  <f>(x,y)  was  constructed  by 
Chen  and  Mei : 


FU)  = If  1/2  (h(V<t>)2  - if2}  dA 

O 


8(<f>  - $_)}  3(4>t,  ~ ) 

+ 1/2  §{h(*R  - ^ da  - §(hcf>a  ^ da 

9.  8* 


- § { h4>  t — } da  + §{h<j>T 
a on  i 


9Ur  - } 


where 


A = region  inside  the  harbor 

§ = line  integral 

4>  = far  field  velocity  potential 

R 
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<J>j  = velocity  potential  of  the  incident  wave 

n = unit  normal  vector  outward  from  region  A 
a = boundary  of  region  A 

<j)  = total  velocity  potential  evaluated  on  boundary  a 

8. 

Proof  was  given  by  Chen  and  Mei  that  the  stationarity  of  this  functional 
is  equivalent  to  the  original  boundary  value  problem. 

9.  The  integral  equation  obtained  from  extremizing  the  functional 
is  solved  by  utilizing  the  finite  element  method.  This  method  is  a 
technique  of  numerical  approximation  which  involves  dividing  a domain 
into  a number  of  nonoverlapping  subdomains  which  are  called  elements. 

10.  The  solution  of  the  problem  is  approximated  within  each  ele- 
ment by  suitable  interpolation  functions  in  terms  of  a finite  number  of 
unknown  parameters.  These  unknown  parameters  are  the  values  of  the 
field  variable  <{>(x,y)  at  a finite  number  of  points  which  are  called 
nodes.  The  relations  for  individual  elements  are  combined  into  a system 
of  equations  for  all  unknown  parameters. 

11.  In  the  region  outside  the  basin,  the  velocity  potentials  are 
solved  analytically  in  terms  of  unknown  coefficients.  The  region  is 
considered  a single  element  with  an  "interpolation  function"  given  by 
Equation  3.  The  infinite  series  is  terminated  at  some  finite  value 
such  that  the  addition  of  further  terms  does  not  significantly  influence 
the  calculated  values  of  <j>(x,y).  The  resulting  equation  is  combined 
with  the  system  of  equations  for  unknown  parameters  at  nodal  points 
within  the  basin  and  this  complete  system  is  solved  using  Gaussian 
elimination  matrix  methods.  n(x,y)  is  related  to  <|>(x,y)  through  the 
linearized  dynamic  free  surface  boundary  condition 


n(x  y)  = _ ± 
mx,y'  g 3t 


The  horizontal  velocity  components  have  the  following  form: 


u(x,y)  = - 


£ 9r>(: 


v(x,y) 


_ £ 3n(x,y) 
w 8y 


MHMIMHIWUMMMBMaM 


r 


12.  The  hybrid  finite  element  method  (so  named  by  Chen  and  Mei 
because  the  method  involves  the  combination  of  analytical  and  finite 
element  numerical  solutions)  is  a steady-state  solution  of  the  boundary 
value  problem.  The  response  of  a harbor  to  an  arbitrary  forcing  func- 
tion can  be  easily  determined  within  the  framework  of  a linearized 
theory. 


PART  III:  RESULTS 


13.  Locations  of  the  three  tanker  terminals  are  shown  in  Figure  1. 
The  finite  element  grid  for  plan  3 is  presented  in  Plate  1.  Plans  IB 
and  3 are  similar  regarding  the  dimensions  of  the  breakwater,  but  rela- 
tive to  the  east  harbor  opening  and  basin  dimensions,  plan  2 is  more 
similar  to  plan  3 than  is  plan  IB  to  plan  3. 

lU.  The  frequency  response  curves  of  the  normalized  maximum  cur- 
rent velocity  (NMCV)  versus  incident  wave  period  at  each  tanker  terminal 
for  plans  IB  and  3 are  shown  in  Plates  2-h.  The  current  velocity  multi- 
plied by  the  incident  wave  amplitude  in  feet  gives  velocity  in  units  of 
feet  per  second.  The  velocities  have  no  vertical  component  or  variation 
since  only  long  waves  are  considered.  Also,  the  maximum  velocity  over 
one  wave  period  is  plotted.  Plates  5-7  are  the  frequency  response  curves 
of  the  wave-height  amplification  factor  (WHAF)  at  each  tanker  terminal. 
WHAF  is  defined  at  a point  inside  the  basin  as  the  wave  height  at  the 
point  divided  by  twice  the  incident  wave  height.  This  definition  of 
amplification  factor  is  traditional  and  is  a result  of  the  fact  that  the 
standing  wave  height  for  a straight  coast  with  no  harbor  would  be  twice 
the  incident  wave  due  to  the  superposition  of  the  incident  and  reflected 
waves. 

15*  The  following  tabulation  gives  peak  values  of  NMCV  and  WHAF 
and  the  resonant  wave  period  for  the  various  Pier  J configurations. 
"Blanks"  in  the  tabulation  reflect  either  the  absence  of  a resonant  con- 
dition or  that  NMCV  or  WHAF  did  not  exceed  1.0.  Several  observations 
regarding  either  the  NMCV  or  WHAF  are  apparent.  For  the  wave  period 
range  under  consideration,  plan  3 exhibits  four  resonant  wave  periods 
(UU,  66,  12U,  and,  nominally,  3^0  sec)  at  each  tanker  terminal.  If  a 
resonant  condition  is  observed  at  a particular  terminal  in  either 
plan  1A,  IB  or  2,  it  is  also  observed  in  plan  3 at  a slightly  different 
resonant  wave  period.  If  the  peak  is  absent  in  plan  3,  it  is  also 
absent  in  the  other  plans.  The  peaks  in  plan  3 are  generally  greater 
than  those  in  the  other  plans  for  the  first  three  resonant  periods. 
However,  for  the  last  resonant  period  (that  covering  a wave  period  range 


from  approximately  2 min  to  well  over  6 min,  and  where  a greater  pro- 
portion of  wave  energy  is  concentrated  per  wave  period  interval) 
plan  3 always  exhibits  a lower  peak  value  than  those  in  the  other  plans. 


STFP  2 (300) 
Plan  1A 


STFP  2 (600) 
Plan  IB 


STFP  3 
Plan  2 


J 

Terminal  1 

i 







1.20  (36) 



3.22  (UU) 

2.08 

(UU) 



— 

— 

— 

2.26  (56) 

— 

3.1*7  (66) 

3.25 

(66) 









1.52  (106) 

1.00  (102) 

2.01  (121*) 

2.93 

(121*) 

2.U5  (280) 

3.93  (285) 

1.93  (260) 

3.22  (268) 

2.50  (350) 

2.75  (31*9) 

1.55  (330) 

2.22 

(31*8) 

Terminal  2 





__ 







2.72  (1*1*) 

1.78 

(1*1*) 





1.65  (59) 



1*. 36  (56) 

1.6l  (56)  ' 

3.1*3  (66) 

2.39 

(66) 









2.55  (106) 

2.55  (106) 

1.91*  (121*) 

3.1*0 

(12U) 

— 

3.1*9  (286) 

— 

3.66  (261*) 

— 

3.21*  (352) 

— 

2.33 

(3U8) 

Terminal  3 

■ • 









1.21*  (1.6) 



3.68  (1*1.) 





1.00  (58) 

1.12  (58) 

1.90  (57) 

3.72  (51*) 

2.79  (66) 

3.80 

(66) 

if  * 









3.1*8  (108) 



1* . 80  (125) 

: ; 

1.78  (280) 

3.69  (281*) 

I.7I*  (260) 

3.32  (270) 

" 

3.12  (360) 

1.72  (330) 

2.25 

(3U8) 

Note:  Numbers  in  parentheses  are  the  wave  period  in  seconds  at  the  peak. 

l6.  Plan  3 evidently  establishes  a harbor  of  sufficient  length, 
width,  and  depth  that  the  first  four  principal  harmonics  can  be  nicely 


generated.  These  harmonics  are  portrayed  in  the  contours  of  the  WHAF 
shown  in  Plates  8-11.  Three  nodal  areas  oriented  parallel  to  the  long 
axis  of  the  basin  are  generated  for  the  first  resonant  wave  period 
(1*4  sec).  The  binodal  and  uninodal  conditions  (each  nodal  area  oriented 
east-west  as  in  the  trinodal  oscillation)  occur  for  the  next  two  reso- 
nant periods  (66  and  124  sec).  Antinodes  form  along  the  inside  of  the 
detached  breakwater  (the  corners  of  which  do  not  particularly  accentuate 
the  oscillation)  and  along  the  face  of  Pier  J.  However,  for  the  6-min 
resonant  period,  a single  antinode  is  formed  in  the  middle  of  the  basin, 
oriented  north-south.  This  is  the  only  harmonic  which  is  similarly 
well  developed  and  oriented  in  all  the  plans. 


t t . 


17.  Plates  12-15  present  the  MMCV  patterns  for  the  resonant  peaks 
of  plan  3.  Comparison  of  these  plates  with  those  from  the  previous 
study*  for  plans  1A,  IB,  and  2 shows  the  currents  are  generally  larger 


- 


in  the  basin  interior  (especially  along  nodal  lines)  for  the  first  three 
resonant  oscillations  (44,  66,  and  124  sec)  of  plan  3 than  for  the  other 
plans.  However,  the  current  patterns  for  plan  3 exhibit  no  particular 
modification  at  the  interior  breakwater  corners  or  high  velocities  at 
the  ends  of  the  breakwater  as  observed  in  the  other  plans.  For  the 
6-min  resonant  condition  of  plan  3,  the  current  pattern  is  similar  to 
that  of  the  other  plans  except  that  as  observed  from  the  preceding 
tabulation,  the  magnitudes  of  the  current  velocities  at  the  tanker  termi- 
nals are  smaller  than  those  in  any  other  plan. 

18.  In  summary,  larger  resonant  oscillations  occur  in  plan  3 
than  in  the  other  plans  for  wave  periods  less  than  2 min.  The  basin 
is  of  sufficient  length,  width,  and  depth  that  three  resonant  periods 
(44,  66,  and  124  sec)  can  be  excited.  The  other  plans  are  of  such  basin 
dimensions  that  one  or  more  of  the  principal  harmonics  are  not  excited 
or  are  excited  to  a lesser  degree  than  that  in  plan  3.  Plan  3,  however, 
excites  the  first  principal  east-west  oscillation  (6  min)  with  lower 
wave  amplifications  and  current  velocities  at  the  tanker  terminals  than 
those  of  any  other  plan. 
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PART  IV:  CONCLUSIONS 
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19-  The  resonant  peaks  of  the  normalized  maximum  current  velocity 
and  wave-height  amplification  factor  do  not  appear  to  preclude  safe 
berthing  at  the  terminal  areas,  given  proper  mooring  practice  and  ade- 
quate mooring  equipment.  The  resonant  response  peaks  for  incident  waves 
with  periods  less  than  2 min  are  generally  larger  for  plan  3 than  for 
the  other  plans.  Plan  3,  however,  exhibits  smaller  response  peaks  for 
the  longer  period  (6  min)  oscillation  than  in  the  other  plans.  Whether 
or  not  incident  waves  in  this  period  range  excite  significant  ship 
motion  depends  upon  ship  characteristics  (length,  width,  and  draft)  and 
mooring  conditions  (type  of  lines,  amount  of  slack,  etc.). 

20.  The  incident  wave  spectrum  and  ship  response  versus  wave 
period  curves  must  be  known  in  addition  to  the  harbor  response  curves 
of  Plates  2-7  to  definitively  determine  whether  or  not  ship  surge 
problems  resulting  from  harbor  resonance  will  occur  in  the  tanker  termi- 
nal basin. 

21.  Until  detailed  information  on  moored  ship  response  as  a func- 
tion of  incident  wave  amplitude  and  frequency  is  obtained,  definitive 
conclusions  on  the  precise  amount  of  ship  motion  cannot  be  made.  How- 
ever, the  oscillations  of  the  tanker  terminal  basin  are  significantly 
smaller  than  the  typical  oscillation  noted  throughout  the  Los  Angeles 
and  Long  Beach  Harbors  complex  in  the  numerical  study  performed  for  the 
entire  harbor  complex.*  Thus,  it  is  reasonable  to  assume,  prior  to 
obtaining  additional  data  on  moored  ship  response,  that  satisfactory 
berthing  conditions  will  obtain  in  the  tanker  terminal  basin. 


* J.  R.  Houston,  "Long  Beach  Harbor  Numerical  Analysis  of  Harbor  Oscil- 
lations; Existing  Conditions  and  Proposed  Improvements,"  Miscellaneous 
Paper  H-76-20,  Report  1,  Sep  1976,  U.  S.  Army  Engineer  Waterways 
Experiment  Station,  CE,  Vicksburg,  Miss. 
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PLATE  1 


NO  LANDFILL 
PLAN  3 
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LEGEND  NOTE  R = WAVE-HEIGHT  AMPLIFICATION  FACTOR 

— NO  LANDFILL  H = WAVE  HEIGHT  , FT  FREQUENCY  RESPONSE 

■—  STFP  2 (600)  Mch=WAVE  HEIGHT  FOR  CLOSED  HARBOR,  FT 

WAVE-HEIGHT  AMPLIFICATION  FACTOR 
62-FT  DEPTH,  TERMINAL  2 


NO  LANDFILL 
62- FT  DEPTH 
44-SEC  WAVE  PERIOD 
SCALE  IN  FEET 


124-SEC  WAVE  PERIOD 


NO  LANDFILL 

62-FT  DEPTH 
44-SEC  WAVE  PERIOD 
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APPENDIX  A:  NOTATION 
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h 
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k 

n 

n 

a 

r 

t 

u 

v 

w 

X 

y 

a 

n 

V 

n 
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Boundary  of  region  A 

Area  of  region  inside  a harbor 

2 

Acceleration  due  to  gravity,  32.2  ft/sec 
Water  depth,  ft 

Hankel  function  of  the  first  kind  of  order  n 
Imaginary  number 
Wave  number,  ft-1 
Integer 

Unit  vector  normal  to  boundary  a 
Spherical  coordinate,  ft 
Time,  sec 

Velocity  in  x-direction,  ft/sec 
Velocity  in  y-direction,  ft/sec 
Angular  velocity,  radians/ sec 
Cartesian  coordinate,  ft 
Cartesian  coordinate,  ft 
Unknown  coefficient 
Gradient,  ft-1 
Wave  amplitude,  ft 

Spherical  coordinate,  degrees 

2 

Total  velocity  potential,  ft  /sec 

2 

Total  velocity  potential  evaluated  on  boundary  a , ft  /sec 

2 

Incident  velocity  potential,  ft  /sec 

2 

Far  field  velocity  potential,  ft  /sec 

2 

Scattered  velocity  potential,  ft  /sec 


In  accordance  with  Bt  70-2-3,  paragraph  6c(l)(b) , 
dated  15  February  1973,  a facsimile  catalog  card 
In  Library  of  Coogreae  format  la  reproduced  below. 


